<p> Draghilev Jacobian/rkfixed Solver </p>

454d0978-e042-4367-9810-3f8fe6c95ef2 11 3 12 decimal

&[DATE] &[TIME] - &[FILENAME]

X0 tmin tmax N D n f rows 1 + : k 1 n range u k el x k el : for b 1 - (f u Jacobian 1 n 1 - n n submatrix * : матрица свободных членов ≡ A f u Jacobian 1 n 1 - 1 n 1 - submatrix : основная(квадратная) матрица ≡ n 2 > det 1 el b A 1 n 1 - 2 n 1 - submatrix augment det - : определителu дополнит.матриц ≡ det n 1 - el A 1 n 1 - 1 n 2 - submatrix b augment det - : det n el A det - : k 2 : k n 2 - ≤ k k 1 + : det k el A 1 n 1 - 1 k 1 - submatrix b A 1 n 1 - k 1 + n 1 - submatrix augment det - : for 41 line det 1 el b : det n el A det : 21 line if dS 1 : i 1 n range D2 i el det i el dS / : for t x D1 D2 : Матрица  правых частей  системы ОДУ ≡ X0 tmin tmax N t x D1 rkfixed 91 line :

x y XY x y / - :

a 816 π * sqrt * : 351.67021

x y f1 1 - 3 / x 2^y 2^+ sqrt * exp 2 y 2^* - x 2^+ a / * - 0.005 + :

Kоординаты начальных точек

f 1 el 1 - 3 / x 1 el (2^x 2 el (2^+ sqrt * exp 2 x 2 el (2^* - x 1 el (2^+ a / * - 0.005 + :

init1 1.78568577383018 021 mat :

init2 1.78568577383018 - 021 mat :

init1 1.78568577382905 021 mat :

init2 1.78568577382905 - 021 mat :

more accurate bounds bound 1.78568577382905 14.24552938284325 21 sys :

Contours are defined by β ... solutions are hyper-reflexive.
'N' defines the fineness ... be generous.
Each contour closes on ≠ length of respective solutions
The best closed contours depend upon 't1max, t2max'. 
Adjust these two manually from the dual observer. tmin 0 : t1max 16149 : N 200 : t2max t1max : 41 sys

B1 init1 tmin t1max N D : B2 init2 tmin t2max N D : 21 sys

The parametrization Human eye by Draghilev method

B1 2 col B1 3 col augment B2 2 col B2 3 col augment 2 1 sys

B1 1 row B1 101 row 21 mat 0 1.786 013 mat 8074.5 14.245 0.039 13 mat 21 mat

dual observers wrt t1max, t2max B1 N row B2 N row 21 mat 16068.255 1.7857 0.0004 - 13 mat 16068.255 1.7857 - 0.0004 13 mat 21 mat

meshed finer tmin 0 : t1max 16149 : N 200 : t2max t1max : 41 sys